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One approach to quantum information processing is to use photons as quantum bits and rely 
on linear optical elements for most operations. However, some optical nonlinearity is necessary 
to enable universal quantum computing [l}{4]. Here, we suggest a circuit-QED approach to linear 
optics quantum computing in the microwave regime, including a deterministic two-photon phase 
gate. Our specific example uses a hybrid quantum system comprising a LC resonator coupled to a 
superconducting flux qubit to implement a nonlinear coupling. Compared to the self-Kerr nonlin- 
earity, we find that our approach has improved tolerance to noise in the qubit while maintaining 
fast operation. 



PACS numbers: 03.67.Lx, 84.40Dc, 85.25Cp 

Linear optics quantum computing (LOQC) has proven 
to be one of the conceptually simplest approaches to 
building novel quantum states and proving the possibility 
of quantum information processing. This approach relies 
on the robustness of linear optical elements, and implic- 
itly requires an optical nonlinearity . Unfortunately, 
progress towards larger scale systems remains challeng- 
ing due to the limits to optical nonlinearities, such as the 
measurement of single photons E] . 

In this letter we suggest recent advances in circuit- 
QED, in which optical and atomic-like systems in the 
microwave domain are explored for their novel quantum 
properties, provides a new paradigm for LOQC [THS]- 
Specifically, using superconducting nonlinearities in the 
form of Josephson junctions and the related quantum de- 
vices such as flux and phase qubits [TTJ [T2] , key elements 
of the LOQC scheme have been realized: the creation 
of microwave photon Fock states [9l [T31 - H5] , controllable 
beam splitters HH] , and single microwave photon de- 
tection [TTJ UH]. In many cases, the photons stored in a 
transmission line-based resonator or inductor-capacitor 
resonator have much better coherence times than the at- 
tached superconducting qubits [TH1 HOJ (32- This sug- 
gests that the main impediment to photon based quan- 
tum computing is the realization of appropriate photon 
nonlinearities to enable two-qubit gates like two-photon 
phase gates, which are sufficient for universal quantum 
computation [T] [5T] . 

The key element for a two-photon phase gate is a two- 
photon nonlinear phase shifter. It imparts a 7r phase on 
any state consisting of two photons, while leaving single 
photon and vacuum states unaffected. In Ref. pQ, this 
scheme is probabilistic and requires additional resources 
to increase the probability of success. Other approaches 
to nonlinear optics quantum computing so far have made 
use of the Kerr nonlinearity [THl [23] . Instead, in this 
Letter we explore the possibility of using a dc SQUID 
[M] to implement a nonlinear coupling between qubit 
and resonator which enables a high fidelity, deterministic 
two-photon nonlinear phase shift in the microwave do- 
main. The nonlinear phase shifter, when combined with 
50/50 beam splitters, can then be used to construct a 




(0 + (d) 




FIG. 1: (a) Implementation of a high-impedance coiled res- 
onator (in blue) coupled to a dc SQUID (red) with an in- 
ductive outer loop. The flux bias lines are in black, (b) 
A simple circuit model of our physical implementation, (c) 
Bottom: Energy levels of the coupled system with a sizeable 
two-photon coupling. Top: The suggested flux bias pulse <j> x 
to implement the nonlinear phase shift; a fast but adiabatic 
sweep and then a very slow variation of the pulse near the 
avoided crossing, (d) Use of two nonlinear phase shifters, 
combined with 50/50 beamsplitters, leads to a deterministic 
two photon phase gate using dual rail logic; storage cavities 
(not shown in (b)) are blue lines. 



two-photon phase gate using so-called dual rail photon 
qubits [HUTU], in which the logical basis {|0) L , cor- 
responds to the existence of a single photon in one of 
two resonator modes (Figure [TJi) . Our approach takes 
best advantage of the relatively long coherence times for 
microwave photons in resonators, and couples only virtu- 
ally to superconducting quantum bit devices, minimizing 
noise and loss due to errors in such devices. When com- 
bined with the aforementioned techniques for Fock state 
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generation and detection, along with dynamically con- 
trolled beam splitters, this provides the final element for 
linear optics-like quantum computing in the microwave 
domain. 

We now outline our approach. We consider photons 
stored in a high-impedance microwave resonator |25j cou- 
pled inductively with strength < \ < 1 to a flux 
superconducting qubit (SQ) in a dc SQUID configu- 
ration (Figure The resonator loops around the 
dc SQUID which results in a nonlinear cosine depen- 
dent interaction between the resonator and qubit. In 
this configuration, we get an effective coupling of the 
form V ~ Ejcos((p + 4>' x ) cos</v, where an external flux 
<f>' x = 2irx&' x /$o is applied to the resonator which conse- 
quently threads the smaller loop of the dc SQUID, $ be- 
ing the superconducting flux quantum. The qubit phase 
variable and the resonator flux are denoted by <p and 
6 r = 27r<Iv/$o respectively. For <p' x ~ tt/2, we see im- 
mediately a nonlinear coupling between the qubit and 
resonator: V ~ Ej4>4> 2 , where two resonator photons can 
be annihilated to produce one qubit excitation, analogous 
to parametric up conversion in \^ systems. This causes 
the two-photon state of the resonator to couple to the 
first excited state of the SQ with strength g 2 (Figure [lj:). 
In essence, in this region, the two-photon state becomes 
slightly qubit-like and acquires some nonlinearity. How- 
ever, the single-photon state, inspite of its coupling to 
the first excitation of the SQ with strength gi, remains 
mostly photon-like because it is far detuned from this 
qubit excitation. At the end of the procedure, this leads 
to an additional phase for the two-photon initial state. 

The noise in the SQ may slightly limit our nonlinear 
phase shift operation. Suppose the first excited state of 
the SQ has a decay rate 7. Although the overall system 
will mostly be in the photon-like regime, there will be a 
non-zero probability for it to decay due to its coupling 
to the lossy qubit. In the limit where the two-photon 
and qubit detuning 8 3> <?2, the two-photon nonlinearity 
goes like g\j8. In the same limit, the two-photon state 
decays approximately at a rate 7<? 2 /<5 2 . Thus, our pro- 
tocol improves as "f/8 becomes small, provided we allow 
gi to become close to 52 in this limit, which is possible 
by making use of <j>' x . In contrast, a Kerr nonlinearity 
scales like g\ /8 3 and the noise scales like -fg 2 / 8 2 , leading 
to more loss at large detuning. 

We now examine a detailed model to support these 
qualitative arguments. In our case, the second resonator 
is not coupled to a SQ and is not shown; we focus on the 
dynamics of the first resonator, which is coupled. The 
quantum Hamiltonian of the system is H = T + V |26j . 
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In addition to <f>' x , an external flux bias <j) x — 27r$ x /$o 
is applied to the outer inductor loop attached to the 
squid. The canonical coordinates of the qubit satisfy 
[<fi,N] — i, where N — q(2e)^ 1 is the number of Cooper 
pairs in the junctions. The operators $ r and q r rep- 
resent quantum fluctuations in flux and charge of the 
resonator satisfying [$ r ,g r ] = iK, and x 1S the fraction 
of the flux $ r threading the squid loop. This inductive 
coupling causes the effective capacitances of the resonator 
and qubit to be modified, and we denote these modified 
values by C r and Cj respectively. Ej is the Josephson 
energy of each junction, while El = $0/ '(4tt 2 £i) repre- 
sents the inductive energy of the qubit due to the bigger 
loop. We can also define an effective charging energy 
of the junction to be Ec — (2e) 2 C7 . We introduce 
another dimensionless parameter fi = 27r"I> " 1 $°, where 
= \J L r uih/2 is the width of quantum fluctuations in 
the resonator flux. In terms of the quantum of conduc- 
tance Go = 2e 2 /h and the characteristic impedance of the 
resonator Z = (Lr/Cr) 1 ^ 2 , we can write /i = \/2ttGqZ. 
Since /x is much smaller than unity, we can expand V in 
powers of x4*r °^ fJ-- Performing the expansion to second 
order, we get H = H r + H q + Vj with 



H„ = 



2Cj 



<T» 2 

— 2Ej cos 



(3) 

















cos 




2 







Vj = xEj[$ r sm( 



X4>1 

2 



cos 



XIQr 

2C T 



corresponding to the resonator, qubit, and interaction 
terms respectively. We remark that asymmetry in the 
two Josephson junctions leads to additional terms, but 
our general linearization approach described below re- 
mains valid, and provides qualitatively similar results. 

The equations (|3| can be used to solve this system nu- 
merically. However, in the regime where El 3> Ej we 
can get some analytical results. First, we linearize the po- 
tential in ([2]) around the classical values of the resonator 
reduced flux <f> c i and the qubit phase (3 c i = — (j>x + f, with 
quantum fluctuations (p r and tp around them. Any non- 
linearity can then be treated as a perturbation. Note the 
following functions of 4> x and <j)' x . 
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E L + Ej[cos<j) x + cos((j) x -<f)' x - X<t>d)] ' 
r = sin/3 c; ; s = sin [/3 ci + <ft f x + x<t>d] , (6) 
t = cos /3 ci ; u = cos [j3 c i +c/>' x + x^ci] ■ (7) 

With the effective inductance of resonator L^ 1 = L~ x + 
(2it/3> ) 2 Ejx 2 u and S r — <& / (2ir)<p r , the resonator and 
qubit Hamiltonians can now be written as 
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with respective frequencies u> — (L r C r ) 1 / 2 and uj q — 
\Ju>c [uil + tuj(t + u)]. We see immediately that chang- 
ing the external fluxes changes these frequencies, and 



hence, the qubit-resonator detuning A 



u>. Intro- 



ducing creation and annihilation operators for the res- 
onator and qubit satisfying [a, a 1 '] = 1 = [b , w] with 

»- v^ (S+t,)! *=- < v'S (S - t,) ' (9) 

the resonator and qubit Hamiltonians become H r = loo) a 
and H q = u) q wb. When the qubit is not linearized as in 
([3]), the potential energy terms can be written as 



V\ = 7]i(a + fit) sin(0 + <f/ x ), 
V 2 = r ]2 (a + a J< ) 2 cos{ct> + (j)' x ), 
V 3 = ir) 3 (a-tf)N, 



(10) 
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with coupling coefficients 771 = 772 = rji/(2Ej), 

and 773 = (r/ihuj)/ (2Ej). The potential that is of rele- 
vance is V2 from which the nonlinear coupling term g 2 is 
seen to be 
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where the matrix element is between the ground and first 
excited qubit states. The size of g 2 is important for the 
success of the nonlinear phase shift protocol. For a given 
value of 771, the Josephson energy Ej of each junction 
cannot be made too large as this will suppress 772 and 
g 2 . Hence, it is desirable to operate the qubit in the flux 
regime where Ej < IOEq ■ The linear coupling coefficient 
771 depends on the characteristic impedance Z of the LC 
circuit implicit in the parameter \i. Therefore, we have 
to implement a resonator with high Z to make g 2 larger. 

After linearization, the flux dependent linear Hamilto- 
nian of the system can be written as 

y 

H L = H r + H q - T^j-QQr + xEju r (p. (14) 

We neglect all higher order nonlinear terms and only con- 
sider the perturbative x^ 2 ) type nonlinearity given by, 
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We can make the rotating wave approximation and write 
Hl in terms of creation and annihilation operators as 

H L = bja ] a + UqPb + gi(ati +a f fe), (16) 

where the linear coupling g\ is given by 



To diagonalize Hl we define new operators c and d as 



/iic + v\d\ b = /j 2 c + v 2 d, 



(18) 



such that [c, &\ = 1 = [d, S] and [c, &} = = [c, d). This 
requires the conditions 

1 = M 2 + M 2 ;mi/4 + ^2=0. (19) 
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The parametrization fix = c 
v 2 = cos 9 satisfies the constraints rtT9|) . Substituting the 
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relations ( 18 1 into Hl and setting the diagonal terms to 
zero, we get a normal mode Hamiltonian Hn — c + 

0,2$ d with energies fii >2 = cj + A/2 (l =F \/l + 45 2 /A 2 ). 

We assume the detuning A > 0. The bare basis states of 
the resonator-qubit system will be denoted by |n) (g> \ q) = 
|n g) , where the first and second labels refer to the quan- 
tum numbers of the resonator and qubit respectively. 
The relevant eigenstates of the Hamiltonian in the new 
basis are number excitations of the operators &c and 
(Pd. Denoting these kets as \CD^, we can write down 
three important eigenstates with energies 0\, 2 , and 
2Qi. They are 

1 10) = cos 6> 1 10) +sin0|Ol) , (20) 

|01) = -sin^ 1 10) +cos0|Ol) , 

|20) = cos 2 6»|20) + V2 cos 6 sin0 |11) + sin 2 6\02) . 

The parameter 9 satisfies tan 29 = —2gi A -1 . For A 
gi, 1 16) -> 1 10), |01) -> |01), |20) -> |20), fi x -> w, and 

2 -> LOq. 

The nonlinearity in our model couples the states 1 20) 
and 1 01) leading to a sizeable avoided crossing in Fig- 
ure^ between the two-photon and qubit levels. In terms 
of the normal mode operators 

V 2 = ?7 2 [cos 2 9 sin 9c n c - cos 3 9c i2 d - 2 cos 9 sin 2 6$cd} 
+2 cos 2 9 sin OtfSd + sin 3 9cA n — sin 2 9 cos 9$ 2 d] + HC, 

where rj 2 = r\ 2 sj\[2. Then the overall Hamiltonian of 
interest becomes H — 0\&c + Vt^d'd + V%. 

Working in the truncated subspace spanned by the 
states {|0) = |06),|o) = |10),|6) = |20),|c) = |01)}, 
we write the Hamiltonian as H = H + V where H = 
O x \a) (a\ + 20i \b) (b\ + 2 \c) (c| and the coupling V = 
Aifla) (6| + 1 6) (o|) + A 2 (|6) (c| + |c) (b\). The parameters 
Ai = ^/2i]' 2 cos 2 9sin9 = r\r( 2 and A 2 = — ^-q^cos^ 9 = 
r 2 ?7 2 . We can adiabatically eliminate the state \a) to find 
an effective Hamiltonian 
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n 2 \c) (c\+rW 2 (\b) (c\ + \c) (b\). 



(21) 



We can use this Hamiltonian to calculate the two- 
photon nonlinearity AT;. For \S'\ = \Q% — 2f2x| ^ l^a^li 
we have 
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where we have associated the nonlinear coupling gi with 
rj^Ti- The nonlinear phase-shift protocol requires initial- 
izing the system in the states 1 10} s» 1 10} and 1 20} ss 1 20} 
with errors that go like g\j A 2 . Then the external fluxes 
cj) x and 4>' x are varied adiabatically so that the state 1 20} 
becomes slightly qubit-like, mostly because of 1 1 1} . After 
accumulating the desired phase, the process is reversed to 
retrieve the photons. For some integer n, we require for a 
total time r g , f s Ni(t)dt = (2n + l)7r. The final outcome 
is th en ^(|00) + |10} + |20}) ^(|00) + 1 10} - |2G». 

In addition to our analytical model, we also diagonal- 
ize the Hamiltonian of the system numerically by work- 
ing in the tensor product space H = H r <S> H„oi the 
resonator and qubit using the Hamiltonian in pi). The 
basis states in the resonator space are the number ex- 
citations \n). The qubit space is written in the basis 
of qubit wavefunctions ip q (<j)) — (4>\q)- We let ft = 1 
and choose u>c = 1 GHz, ujj = 5 GHz, lul = 3w,/, 
and to = 2.225 GHz. The characteristic impedance 
Z w 449 fl. We choose a \ = 0.17, representing an 
easily achievable mutual inductance, from which follow 
77i = 400 MHz, 772 = 16 MHz, and 773 = 89 MHz. 



6 (GHz) 




FIG. 2: (a) Schematic of the system bare energy levels and 
couplings, (b) Contour plot of detuning 5 and | c/2 1 with the 
on and off points marked in green and red. The on point is 
chosen such that the (72 is maximized, (c) Top: The coupling 
gi/W and (72, with the on and off fluxes shown. Bottom: 
The bare frequencies 2uj and ui q obtained from the analyti- 
cal model (dashed) and numerical results (solid). The over- 
all qubit-resonator interaction leads to a splitting of approxi- 
mately 10 MHz. 

Within the linear model, we can take into account loss 
of the qubit excited state at a rate 7. The expectation 
value of the operator &c can then be calculated in the 
dressed eigenstates 1 10} and 1 20} . In the limit A> \g\\, 
it will be found that the two photon sector decays at a 
rate jg 2 /A 2 = "fgl/(5 + lu) 2 . For \5\ 3> I52I, there is also 



a decay rate jg% / S 2 due to the nonlinear coupling of the 
two-photon sector with the qubit. Including the cavity 
decay rate k, the total decay rate of the photons becomes 
r = K + 1 g 2 /A 2 +1 g 2 /5 2 . 
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FIG. 3: (a) A plot of the dimensionless dynamic loss Ld 
for k — \ KHz, 7 = 100k and e 2 = 0.01. The detuning 
-536 MHz < 8m < -41 MHz. (b) The static loss L s in pur- 
ple, and the static loss without the effect of the cavity decay 
rate k in green. 

Assuming g 2 is time independent for simplicity, adia- 
baticity requires g%\S\ 2 (5 2 +4g% )~ 3 <C 1 . We can set this 
equal to some e 2 <C 1 and solve for 
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— wdS, (23) 

St (<5 2 +4 5 2 )s 



which is the time taken to go from \8i\ ^> |<?2| at t = to 
smaller values of detuning with a minimum 8 m . Taking 
6 <72 adiabatically will take time -77 ~ (7|^2 1 ) 1 so 
that the total time of the protocol r g > 2t^. The total 
dynamic loss during the process is given by 



L d = 



Si 



(S + uj) 2 



1.92 



3 dS. (24) 



When the detuning is held at 6 m for a time t s , the static 
loss L s is approximately 
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(25) 



It can be seen from the Figure (3]d that the photon loss 
for large detuning is only limited by the cavity decay rate 
k. We require the total loss L = L4 + L s <C 1. 

The protocol might also be limited by dephasing of the 
qubit due to flux noise [2"TH2l)] . The average slopes of the 
single and two-photon energy levels with respect to the 
reduced flux (j) x are approximately 50 MHz and 100 MHz 
respectively, while the slope of the qubit energy level is 
at most 1 GHz for the parameters chosen. However, the 
exact loss due to dephasing depends on the flux noise 
amplitude [30 | [31 ] . 

In conclusion, we have demonstrated that by appro- 
priately tuning the two control fluxes, the nonlinear cou- 
pling enables a two-photon nonlinear phase shift opera- 
tion with loss at large detuning limited only by the cavity 
quality factor. This is highly desirable compared to the 
self-Kerr nonlinearity which leads to more loss at large 
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detunings. Furthermore, our approach may be adaptable 
to recent ultra-high quality factor resonators enabling 
nonlinear optics quantum computing in a fully engineered 
system [52] . 
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